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Abstract 

A new multifermion formulation of lattice QCD is proposed. The model is 



o 

m 

> 
o 

■ free of spectrum doubling and preserves all nonanomalous chiral symmetries 

, up to exponentially small corrections. It is argued that a small number of 

QQ I fermion fields may provide a good approximation making computer simula- 

0> ' tions feasible. 

1 Introduction 

It was believed for a long time that any local regularization of a theory which 
posesses a chiral symmetry breaks a manifest chiral invariance. In the context of 
^ ' lattice models this statement is the essence of the Nielsen-Ninomiya "no-go" theorem 

|T[]. However in our paper it was shown that anomaly free models, such as the 
Standard Model, do allow a local chiral invariant regularization if one introduces 
the regularized action containing infinite number of fermion fields. Application of 
this idea to lattice theories lead to the models which have no spectrum doubling and 
posess exact chiral invariance in the continuum limit 0, Q]. The method was also 
checked by nonperturbative simulations of some 2d models 0. R.Narayanan and 
H.Neuberger showed that this approach is closely related to the five dimensional 
domain wall model proposed by D.Kaplan 0. Introduction of an infinite series of 
fermion fields may be considered as a discretization of the fifth dimension. 

Although originally it might seem that considering an infinite series of fermions 
is too high price, by now it is widely accepted that it is the only way to preserve 
the exact chiral symmetry without tuning additional parameters. Many present 
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approaches like overlapp formalism or different versions of domain wall fermions 
§, dig use this idea. 

The main question is whether these approaches allow efficient nonperturbative 
calculations. Considering realistic chiral gauge models on the lattice seems at present 
to be too hard for computer simulations. On the top of the usual problem of simu- 
lating fermion determinants there is an additional difficulty, as a chiral determinant 
being complex does not allow straightforward Monte-Carlo simulations. So it seems 
reasonable to start with the study of vectorial models like QCD and to try to deal 
with the breaking of the global chiral symmetry. 

A naive introduction of the Wilson term in QCD to remove the degeneracy of 
fermion spectrum breaks global chiral invariance of the action leading to appearance 
of quark mass counterterms as well as order a chirality breaking corrections. Apart 
from esthetical objections it raises practical questions for numerical simulations. 
Whereas the Abelian chiral symmetry is indeed broken by the triangle anomaly, the 
SU{2) chiral invariance is not affected by anomalies and is believed to be broken 
spontaneously. Explicit breaking of the SU (2) symmetry by the Wilson term leads 
to a nonzero mass of 7r-meson which is supposed to be a pceudogoldstone boson. 

To supress these effects one may try to modify the Wilson action by adding the 



terms killing the order a corrections [11]. It partially solves the problem, but the 



modified action makes computer simulations more difficult. Extensive study of this 
approach was carried out in the papers (see [jl2| and references within). 

It seems preferrable to use some approach which preserves, at least approxi- 
mately, the global chiral symmetry of the theory. One possibility, which is discussed 
now, is to use the actions, satisfying the Ginzparg- Wilson relation . It was shown, 
that if the model satisfies this relation, a number of results which one expects in a 
chiral invariant theory follows , [|l^ . Moreover in this case some modified chiral 



symmetry is present |TB|. Examples of such actions were given in the papers |T7 



||16|| . They are typically nonlinear and formally nonlocal functionals, but nonlocal 
effects are damped exponentially. Some qualitative confirmation of the relevance of 
this approach for low energy QCD was presented in the paper 



The multifermion models like truncated overlapp |TP| or domain wall fermions 
T0[] were also applied to the study of QCD and first numerical simulations ||21| 
gave promising results. The chirality breaking effects, in particular fermion mass 
renormalization, were also discussed |Q, p3| . 

In the present paper I propose a new version of a multifermion formulation 
of QCD. The model has no species doubling, does not require any fine tuning to 
calculate arbitrary gauge invariant amplitudes and does not contain any dimensional 
parameter except for a lattice spacing. For a finite lattice spacing all chirality 
breaking effects are supressed exponentially. Although the model includes an infinite 
number of fermion species, the convergence of the corresponding series is very fast, 
and one can hope that it may be truncated by a small finite number of terms. 

The paper is organized as follows. In the Section 1 we formulate a general idea 
of the method and prove that our effective action produces the fermion determinant 
differing from the gauge invariant QCD determinant by exponentially small terms. 
In the second section we study the chiral property of the model and demonstrate 
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the absence of a perturbative quark mass renormalization as well as exponential 
supression of all cliirality breaking terms. In the Section 3 it is shown that all 
nonanomalous chiral currents are conserved whereas the usual Abelian anomaly is 
present. In Section 4 we discuss possible applications of the model to numerical 
simulations. 



2 Representation of the QCD quark determinant 
in the multifermion model. 

To give an idea of the method we start with the continuum model describing a gauge 
invariant interaction of two generations of spinor fields Xn, n = 1 . . . oo with the 
masses proportional to the number n, rUn = mn. The following equality holds 



exp{- EE/ xt{x)iD + mn)Xn{x)dx}dxtdXn= (1) 

fc=l,2n=l 

/oo ^ 
exp{- E / Xn{x){~D'^ + m'^v?)xn{x)dx]dxndXn = 

n=l 

cncoth(^)A 

Here Xn are spinor fields with the Grassmanian parity (—1)", -D is a QCD Dirac 
operator and Di are it's eigenvalues. C is a field independent constant. 
To derive this equation we note that the integral over x is equal to 



I =f[ det{-D^ + 



2^2^-1) 



m n 



n+l 



n=l 



This equation may be rewritten as follows 

oo 

I = U^Mi: ln(A' + mV)(-l)"+i} (3) 



n=l 



Summation in the exponent @ may be done explicitely, giving the result 

oo TT I D I 

y \n{D^ + mV)(-l)"+^ = ln(coth + In \DA (4) 

^1 2m 

This proves the eq.( |l|). 

It follows from the eq.(|l]) that up to a nonessential constant factor 



m— >0 

and 



lim^J = det(D) (5) 
lim / = 1 (6) 



m— >oo 



3 



These equations strongly suggest that if we apply analogous procedure to lattice 
QCD, and instead of the mass term introduce the Wilson term W, we get a theory 
which for small values of W produces the QCD massless quark determinant and for 
large values of W (doublers region) is trivial. As follows from eq.(|^) the doublers 
contribution is supressed by the factor exp{ — |^}, so we expect that all chirality 
breaking effects are also supressed exponentially. Below we justify this conjecture. 

Note that in distinction of our previous approach 0, ^ as well as different 
versions of domain wall fermions, which include heavy states in the domain p ~ 0, 
in the present version only massless states are present in this region, and heavy 
states are confined to the doublers region p ~ vra"^. No additional mass parameters 
appear in the theory. 

Our claim is that the QCD determinant may be presented as the path integral 
of the exponent of the following action: 



oo 



- ^ ^ -rix)^b,iD, + D;)-nKaD;D,]M^) (7) 

^iM n=— 00,717^0 

D^ = ^[U^{x)ij{x + a^)-ij{x)] (8) 

Here the first term represents the chiral invariant part of the action and the second 
term is the Wilson term multiplied by n. In this equation k is a dimensionless 
parameter, which we choose in the interval < k < 1. The fields ip"' have a 
Grassmanian parity (—1)"". 

In distinction of the eq.(|l|) where we considered two generations of xt fields 
(fc = 1, 2), and 1 < n < oo, in eq.( |^ we introduced one generation, but allow n to 
change in the interval — oo<'ri<oo, n ^ 0. For a zero bare quark mass it makes 
no difference, but , as it was pointed out to me by A.Ukawa, for nonzero quark mass 
the choice —oo<n<oo, n is more convinient. 

Let us consider a fermion loop diagram with L vertices, which may be presented 
by the following integral 

ndqi . . .gL-i) = / d'p (-l)"Tr [G(p)r,(p, gi)G(p+gi)r,(p, gs) • • •] = 

(9) 

-1 oo 



2 [ y( i^n Foip,q) + n'F,{p,q) + 



Here Fj denotes the interaction vertices and qi is a total momentum entering the 
corresponding vertex. The fermion propagators G[j)) look as follows 

^ , ^ is + mn , , 

= 2^ 2 2 10 
gZ _|_ jYi^^^z 



a^^sm[{p + q^ + ... + q^)^a] (11) 
^^(l-cos(p + gi + ... + g%) (12) 
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The eq.(^) may be rewritten in the following form 

A,(p,q) 



n, = 2/"\d>x:E(-iryS (13) 

J —na ^ , n. 1 oi nil II 



,^0 n=l + 

Here Ai are some functions of p, g, which do not depend on n. They satisfy the 
equations 

n sl = Fo{p.q) (14) 

I i=0,ijtl 

I 

The summation over n can be done as before, giving the result 

n, = r' rfvi- 1:' — + ^] (15) 

^ '-0 mi J sf sinh{7r J nil ^' 



One sees that for small p+qi + . . . + qi the first term is damped exponentially and the 
second one describes the amplitude for massless quarks. For p + qi + . . . + qi ~ vra"^ 
the first term is 

A{p,q) 

(16) 

and exactly compensates the second one providing the doublers supression. 

Now we shall study eqs. (H- |15D in more details and prove that the amplitudes (H) 
coincide with the gauge invariant QCD amplitudes up to the terms which decrease 
exponentially when a — 0. 

Firstly we consider the case when all momenta qi are external. Then we can 
assume that 

\qi\<ea-\ e<l (17) 
The integration domain in the eq.(|^) can be separated into three parts: 

Vi : \p\ < ea-^; V2 : Lea"^ < \p\ < (vr - Le)a"^; (18) 

V3 : (vr - Le)a~^ < \p\ < TTa'^ 
In the region V2 the function 

- {J2 sin(P/.a)(l + cot(p^a)(gi + . . . + qi)a)))~^ (19) 

may be expanded in terms of qia. It is easy to see that the the functions Ai also may 
be expanded in terms of qia. Therefore integrating the second term in the eq. (|T5|) 
over p one gets a local polynomial in qi. Analogous reasonings show that the first 
term in the eq.(|l5D also contributes a local polynomial in g;. 

In the domain V3 we cannot use the expansion (|T^) as in this region sin(p^a) ~ 0. 



However the complete r.h.s. of the eq.([r5|) still allows the expansion in terms of qia. 



5 



Indeed, in this region sf < mf, and the r.h.s. of the eq(^) may be presented as 
follows 

The terms ~ s^"^ cancel and the remaining terms may be expanded in terms of qia. 

In the region Vi we use the representation ([T5|) for the amplitude (|]). In this 
region \a{p — qi)\ < 2e and therefore 

J2 sin^[a(p - qi)f,] >[J2l- cos(p - qi)f,a]'^ (21) 



Hence the first term in the eq.(|T5D is supressed exponentially by the factor 



sinh(7ry S;^m; ^) ^'^ 

The second term gives 



^ exp{--} (22) 



wt^-r:^'^^ (23) 



where the eq.(14) was used. This is exactly the expression which one would get 
in the chiral invariant theory with no Wilson term. In the limit a — this term 
reproduces the diagrams of the continuum QCD. It is worthwhile to notice that if 
only this term were present, the gauge invariance could be broken, as cutting the 
Brillouin zone is not a gauge invariant procedure. In our construction the gauge 
invariance is preserved at any stage of calculations. Possible noninvariant terms 
arising from cutting the Brillouin zone in the eq. ([25|) are compensated by local 
counterterms generated by the integral over V2, V3. It would be wrong just to drop 
the integration over V2, V3. The local counterterms coming from these integrals are 
crucial for maintaining the gauge invariance of the amplitudes. So we proved that 
for one loop diagrams, when all qi are external momenta, qi < ea~^, the amplitude 
(^ differs from the gauge invariant QCD amplitude by exponentially small terms. 
The value of the correction term for the fixed a depends on external momenta q. 
In particular for g = it is zero for any a. The rate of the exponential damping 
depends also on the parameter k, and one may have a temptation to take k very 
close to zero. It is clear however that if one puts k = 0, the correction term indeed 
disappear, but the doubler states reappear. To determine possible values of k we 
note that considering the contribution of the region V3 we used the expansion in 
terms of qa. It is easy to see that for p ~ vra"^ the effective expansion parameter 
is — . So the expansion to make sense we need k, > qa. In particular we may take 
K = e. 

In the next section we consider the multiloop diagrams and prove that they also 
reproduce chiral invariant QCD amplitudes up to exponentially small corrections. 
To do that we have to study more closely the global chiral invariance of our model. 
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3 Global chiral invariance of the model. 



In this section we study the invariance of our model under global chiral transforma- 
tions. We show, that as expected in a chiral invariant theory, a perturbative quark 
mass renormalization is absent. All chirality breaking effects are supressed exponen- 
tially. The divergency of any amplitude including an SU{2) chiral current vertex, 
or any other chiral current, corresponding to an anomaly free group is zero up to 
exponentially small corrections. At the same time the divergency of the Abelian 
axial current is not zero, producing the usual chiral anomaly. 

To give an idea of the mechanism which provides the supression of chirality 
breaking effects we firstly consider a simplified case. All chirality breaking effects 
arise due to the presence of the Wilson term and therefore vanish when n = 0. 
Let us suppose that some chirality breaking term, for example the fermion mass 
renormalization, is simply proportional to n: 5mn = g'^M{g)n. 

Introducing such counterterm is equvalent to shifting the Wilson term by g^Mn. 
To get the expression for the QCD determinant we still may use the eq.s( ^ |15D, 
changing everywhere m/ by mi + g'^M. It is sufficient to consider the case when 
\qi\ < ea~^ and the integration domain in the eq.( [1^) is restricted to (— ea~^, ea~^), 
as it was shown above that larger values oip produce only local counterterms needed 
to maintain the gauge invariance. The expression for the QCD amplitude acquires 
the form: 



,4 iMP^Q) 



= , dM^^ - (24) 



L-l 

-E 

/=o {mi + g'^M)^sfsinh{7r^sf{mi + g^M)-^y 
The sum of the first terms in the eq. (|2^) due to eq.(|l4D is equal to 

">V#^^ (25) 

which coincides with the result obtained before for the model without mass countert- 
erm insertions. Expanding the second term into Taylor series over g one sees that all 
the coefficients are exponentially damped, being proportional to exp{— Try^m^^^} < 
exp{— 7r(Ke)^^}. So for the simplified case we proved our statement. 

However in general chirality breaking terms may have more complicated depen- 
dence on n. Below we shall show that our statement remains true in a general case 
as well, although the proof is more involved. 

To prove the chiral invariance of the determinant it is sufficient to show that 
it's derivative with respect to k, is proportional to some chiral invariant amplitude. 
Note that this derivative is not necessarily zero as chiral invariant counterterms may 
also depend on k. We demonstrate that it is indeed true up to exponentially small 
terms. 

We study the following object 

/ = ^( Y: (-irTrlnDO= (26) 

n=—oo,ny^Q 
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= Y.i-'^T Tr ( nW{D + nnW)-^) 

n 

Here W is the usual Wilson term. This equation may be rewritten as follows 

^ = E(-l)"Tr(niyG„) (27) 

n 

wher Gn is the Green function of the Wilson-Dirac operator in the external gauge 
field. This function satisfies the chiral Ward identities which in the continuum 
notations may be written as follows 

Gnix, t/)75 + 'j^Gnix, y) = 2inK J Gn{x, z)-i^W{z)Gn{z, y)dz (28) 

75^^75 = -G-n (29) 

The Green function G„ may be presented as a sum Gn = + where G^ 
is a chiral symmetric part which anticommutes with 75 and is even in n. is a 
chirality breaking part which commutes with 75 and is odd in n. It follows from 
eq(^ that 

Gn{x,y) = inK'y5 / Gn{x, z)'y5W{z)Gn{z,y)dz (30) 



Using these properties of G„ one can rewrite the eq(^) in the form 

00 „ 

/= Yl {-iTn^ii^ Tr{W{x)-f5Gn{x,z)-f5W{z)Gn{z,x))dxdz (31) 

n=— oo,n^0 

The eq(^) is written for one loop diagrams. To get the multiloop diagrams one has 
to integrate it with the exponent of the gluon action. In terms of Feynman diagrams 
it is represented by the fermion loops with two insertions of the Wilson term W. 

We shall start with the diagrams including one fermion loop and arbitrary num- 
ber of gluon vertices and lines. 

The analytic expression corresponding to the eq(PTD for such diagrams may be 
written as follows 

._.xn 2 (-^0 + n^Ki + . . .)F{qr . . . qL-i)dpdqrdqL^i 
nMnJ {sl + mW)...{sU+ml_,n-) ^ 

Here we wrote the expression for the diagram with L quark- gluon vertices. The 
notations are as follows: 



*o = V(II^(P))' "^o = "^(p), (33) 



= v(II^J(P + 9i + • • ■ + qi)), mi = m{p + +qi + . . . + qi) 



The momenta qi,...qr are external and satisfy the condition \qia\ < e. The mo- 
menta qr+i, ■ ■ ■ qi-i correspond to virtual gluons and may acquire arbitrary values 
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in the interval (— ^, The function F{qr, . . . Ql-i) depends only on virtual gluon 
momenta corresponding to gluon propagators and vertices. 
We use the following decomposition of the integrand 

Ko + n'K^ + ... ^ Bi 



E (35) 



Here Bi and are n-independent functions satisfying the relations 

E^^n^H^o, ... E^^n^^Ho (36) 

E^.n^? = i> ••• E^^n^^=o (37) 

k ij^k k ij^k 

We assume that all momenta si are different. It is not true for the diagrams with 
self energy insertions, but it will be clear from the final result that our conclusions 
are valid for this case as well. 

Using these decompositions we may rewrite the integrand in eq(^) in the form 

_ oo R 4 

V (-l)V(V . o J ( E 2 2 . )= (38) 



' ^2^2 ^2^2 V ^2 I ^2^2 ^2 , ^2^2- 



Performing the summation over n one gets 

^ s%mf - sfm% siuii smh(7rs/m^ ) s^mfc smh^vrsfcrn^ ) 

All the momenta s;, > / > r are small and therefore the first term gives exponen- 
tially small corrections. If some of the momenta are also small, the corresponding 
term is also supressed exponentially. So it is sufficient to consider the eq (|39|) when 
all Sk are big: e < |a(p + . . .qk)\ < vr. 

The remaining terms may be expanded in terms of mfsj'^. (We recall that qia < t 
and therefore mi <^ si. 

h= E -u E^(i + ^ + ---) (40) 

k=r+i Skml smh(7rsfcmfc ) sf sfmi 
Let us consider the sum over /. Due to eq(^) the first term is equal to 

i:§ = n=-'?«..-L./2 (41) 

1=0 -bQ. . . 
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where 

_ Tr(gorigir2 . . .) 
ll«=o 

is the manifestly chiral invariant amphtude corresponding to the theory without the 
Wilson term. 

The first factor in the eq(^Dp may be expanded in terms of aqi where qi, I = 
1, . . . r are "small" external momenta. Being integrated over p, q^+i . . . qi-i the 
first term in the eq(^) is equal to the manifestly chiral invariant amplitude 
integrated with the chiral invariant function and multiplied by a polynomial of 
external momenta. Such a term is obviously chiral invariant. 

To analyse the next terms we use the relations which follow from the formal 



expansion of the eq (plf) in terms of n 



E§(l-"4 + ---)= (43) 

1=0 



1=0 1=0 1=0 

Using these relations we can rewrite the second term in the r.h.s. of the eq (|40D in 
the form 

k=r+i smni^TTSfcm;^ ) 



The first term in brackets is equal to 

-^1 nfc=r+l ^1 tO 



(45) 



The factor KiKf^^ Uk si may be expanded in a power series in aqi, I = 1, . . . r. 
Coefficients in this series are not singular at p = 0. Therefore being integrated 
over p,qk, k = r + l,...L — 1 this term is equal to a manifestly chiral invariant 
amplitude Jq multiplied by a polynomial in qi. In the second term in the eq(^4D we 
can expand mfs^'^ in terms of aqi, I = 1, . . .r. Again it is easy to see that the 
coefficients of this expansion are not singular at p = 0. Hence being integrated over 
p, qk this term produces only local counterterms. 

The next order terms in the eq(^) are analyzed in the same way and lead to 
the same conclusion. It proves that the derivative over k of the diagrams under 
consideration is equal to a chiral invariant amplitude plus a local gauge invariant 
polynomial up to exponentially small corrections. Hence the QCD determinant is 
chiral invariant. 

These reasonings are easily extended to the diagrams including more than one 
fermion loop. A diagram with two fermion loops may be presented as a product of 
two diagrams including one fermion loop integrated with some function depending 
on the virtual gluon momenta. For each one loop diagram one has to use the 
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decompositions ( p4| , [36| ). Multiplying the corresponding expressions and using the 
representation 



{sf + mfn'^){sl + mln'^) slmf — mfsl) sf + mfn'^ s\ + m\in? 

one reduces the problem to the one considered above. Diagrams with arbitrary 
number of fermion loops are analyzed in the same way. 

It completes the proof of the chiral invariance of any gauge invariant amplitude 
in our model. 



4 Chiral currents conservation. Anomalies. 

In this section we prove that the currents assosiated with nonanomalous global chiral 
symmetries are conserved up to exponentially small terms. 

Let us consider the action (|^) with the source term for a chiral current: 

/v^ = / + ^y»j;(a;) (47) 
Here is the axial current which in the continuum case is equal to 

oo 

Jl{^)= E M^h.l.T'i^nix) (48) 

n=— oo,n^0 

being the generators of the chiral symmetry group. is the source term. In the 
continuum chiral invariant theory this current is covariantly conserved 

d,Jl + gt^^'jX = (49) 

For the lattice action ( ^7D including the Wilson term the corresponding equation 
is modified. Firstly, the Neuther current has a contribution from the Wilson term, 
and secondly, as the Wilson term breaks the chiral symmetry, the r.h.s. is nonzero. 
The analogue of the eq.(^) for the action (^) may be written as follows 

oo 

d,Jt + gt'^''jX= E nJ2Mx)W{x)Mx) (50) 

n=-oo, n^O X 

where 9 is a lattice derivative. The r.h.s. of eq. (|50|) is due to the Wilson term and 
in the quantum theory is responsible for the chiral anomaly. A particular form of 
W is not essential for our discussion. It is important that the r.h.s. of eq.(|50|) is 
bilinear in ip and proportional to n. 

Being integrated with the exponent of the action (|47|) the eq . (|50D leads to the 
chiral Ward identity broken by the Wilson term 

a-' sm{k^a)T'lXX {k, gi • • • gn) + (51) 
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c?Vp,::r;„(^ - p, gi • • • qn)v^ip) = E 

n=— oo, n^O 

Here Wn stands for the sum of diagrams generated by the r.h.s. of the eg. (|50|) . 

We shall use the differential form of this identity, which is obtained by differen- 
tiating it over and putting k = 0: 

r'tv::X(o,gi...g„)+ (52) 



. 00 
+gt'^'j d'pd,r;--xik -p,q,... gn)U=oV?(p) = E 



Here is the important difference between anomalous and nonanomalous theories. 



Passing from the eq. ([5T|) to the differential identity (E% we assumed that ^ is not 



singular at = 0. It is not true for anomalous models: the anomalous triangle 
diagram is singular at k = 0. Therefore the eq.(^2|) is equivalent to the identity ([51]) 
only if anomaly is absent. 

Now we shall show that the r.h.s. of the eq. ([5^) is exponentially small. To 



prove the absence of anomaly it is sufficient to consider one loop diagrams. The 
diagrams representing the r.h.s. of eq.(^) are obtained from the QCD amplitudes 
which we considered in the previous section by insertions of the additional vertices 
n^^\k=o = nfx, corresponding to zero external momentum. In other words they may 
be considered as momentum dependent counterterms which shift the Wilson term 
nm — > n{m + fi). The sum over n gives for such diagrams the result analogous to 
the eq.dl) 

= r' <i*p|ii%i) - (53) 



^2 



L-1 



- E -Mp^ci) J ^ 

1=0 {mi + fi)Jsfsmh{'iT J 3^(1711 + fi) ^) 



where c.t. stands for local counterterms. As was explained above the sum of the 
first term and local counterterms is equal to the gauge invariant QCD amplitude 
without insertions of fi and the remaining terms are proportional to exp{ — j^^^}. 
The diagrams with k insertions of fi are given by the k-th term of the formal Taylor 
expansion of this exponent in terms of /i near the point /i = 0. All these terms 
are exponentially small. So we proved that to any nonanomalous global chiral 
symmetry corresponds a current which is covariantly conserved up to exponentially 
small corrections. At the same time the global U{1) anomaly is present as it was 
explained before the differential Ward identity (|52D does not hold in this case. 



5 Discussion 

In the previous sections we demonstrated that the lattice QCD determinant may be 
presented by a path integral of a local action, which is free of spectrum degeneracy 
and reproduces gauge invariant amplitudes of massless QCD up to exponentially 
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small corrections. Nonanomalous global chiral symmetries, in particular SU{2) chi- 
ral invariance are also preserved up to exponentially small corrections. No cliirality 
breaking counterterms are needed and the model provides automatically 0(a) im- 
provement. 

Our procedure may be also used in the case when a nonzero bare quark mass is 
present. In this case one has to consider the sums of the form 

oo /I \ n 

E 2 / ^ (54) 

where mo is the bare quark mass. Representing this sum by the integral 

rdz (55) 



L 



c cos{t:z)[si + {miz + mo)^ 

where the contour C encloses the real axis except for a vicinity of the point z — 
one sees that this integral is equal to 

0(exp{— TTSim^"^}) (56) 



Of course in this case the chiral current conservation is broken by the mass term 
and chirality breaking corrections ~ mo are present. 

Om^ effective action includes an infinite series of fields, and for practical simu- 
lations it is important to know how sensitive is the procedure to cutting this series 
by some finite number of terms. Fist of all we note that cutting the series does not 
spoil the gauge invariance of the model. So all the arguments using this invariance 
are still applicable. Due to the fast convergence of the series we expect that a small 
number of terms may provide a good approximation. This number depends on the 
value of external momenta g^, which are relevant for the process under considera- 
tion. As we showed above the integration over momenta bigger than qi produces 
local counterterms. Therefore it is sufficient to study the convergence of the series 
for \p\ < \qi\. The parameter which determines the value of the correction term is 

^J^mJ^ — b. For external momenta \qa\ < 1 this parameter is > 2 (for k — 1). 

The following simple observation may be useful for numerical simulations. For 
large values of b the correction term is small, however to approximate it by a finite 
series one needs to cut it by some number 3> 6. We recall however that in our 
scheme the QCD determinant is presented as a product of the determinants corre- 
sponding to the fields ipn, < n < oo, n ^ 0. To improve the convergence 
for large b one can cut the series asymmetrically — < n < A^-l- 1. Then for b ^ N 
and even A^ the contribution of the fields with negative n is close to zero, whereas 
the contribution of the fields with positive n produces the correct result. 

Numerical estimates of the correction term for go < 1/4, k, — l/A show that 
cutting the series by — 2 < n < 3 gives the relative correction ~ 0.02. Of course it is 
a rather crude and naive estimate and only real simulations may prove the efficiency 
of the method. However it raises a hope that relatively small amount of fields may 
be sufficient to have a strong supression of chirality breaking effects. We recall also 
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that only half of the ipn fields are fermions which is important for simulations with 
dynamical quarks. Having in mind that these simulations do not require a fine 
tuning of mass parameters and introduction of "improving" terms, one may think 
that this model provides a competitive method of simulations in lattice QCD. 
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